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We present evidence about a critical behavior of AD compact QED (CQED) pure gauge theory. Regularizing 
the theory on lattices homotopic to a sphere, we present evidence for a critical, i.e. second order like behavior at 
the deconfinement phase transition for certain values of the coupling parameter 7. 



1. Introduction 

Four dimensional CQED pure gauge theory is a 
prototype model of gauge interactions of abelian 
nature. It has a deconfinement phase transition 
at values of the gauge coupling constant of order 
1, separating confinement and Coulomb phases of 
the theory. The order of the deconfinement phase 
transition in the /3 — 7 plane of couplings 

5= -^[/3cos([/p) + 7cos(2[/p)] (1) 
p 

has been studied for long time. In the early times 
of LGT the transition on the Wilson line (7 = 0) 
was thought to be continuous. Later the tricrit- 
ical point (TCP) was suggested [1] at a slightly 
negative value of 7tc'p = —0.11, predicting that 
the deconfinement phase transition on the Wilson 
line is of weak first order. 

The physically interesting question, whether 
there really exists a critical point in 4_D CQED 
and if so, what are its properties, has been left 
unanswered throughout the years. One can ap- 
proach this question by a determination of the 
correlation length critical exponent i>. A value 
of = i can be expected assuming a mean field 
behavior. A value of i' = jj would correspond to 
a discontinuous phase transition. Other and dif- 
ferent values of i> would be indicative of an "un- 
expected new" behavior. The existence of such 
a topological disorder driven fixed point could 
have implications for the continuum limit of 4_D 
field theories containing compact U(l) gauge de- 



grees of freedom, especially when fermions are in- 
cluded. 

Although a critical transition is expected to ex- 
ist at least at sufficiently large negative 7-values 
[1,2], this was never confirmed in the stochastic 
evaluation of the path integral. In fig. 1 we dis- 
play the probability distribution of the normal- 
ized plaquette expectation value evaluated close 
to the deconfinement phase transition at a 7- value 
of 7 = —0.2 on a S'* lattice with p.b.c. As can be 
seen, double peak structure is present and actu- 
ally can be observed up to quite sizable negative 
values of 7 [1,3]. The appearance of such signals 
may cause doubts about the order of the transi- 
tion and excludes a meaningful finite size scaling 
(ESS) analysis of the critical behavior of CQED. 

We have argued [4], that the disturbing effects 
may be due to an inappropriate choice of the 
boundary conditions of the system. Therefore 
there is need for a improved lattice regularization 
of CQED. Almost all previous numerical simula- 
tions of the CQED path integral have been per- 
formed on a torus with p.b.c. We proposed to use 
lattices topologically homotopic to a sphere. On 
such a lattice any monopole loop, closed on the 
dual lattice, can be contracted to a point, while on 
a torus this cannot be achieved for all loops. The 
actual construction of the lattice was performed 
by considering the manifold given by the surface 
of a 5-dimensional hypercube with linear extent 
L, denoted SH(L) here. Simulations on SH(L) 
on the Wilson line demonstrated somewhat sur- 



prisingly the absence of metastability signals. 
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Figure 1. P(e) on torus at j 



-0.2. 



However, in that study the plaquette terms of 
the lattice action of (1) were discretized in the 
most naive way. Plaquettes in SH(L) were as- 
signed equal weight to the action function. Thus, 
although the global topological properties of a 
sphere were implemented, effects related to the 
curvature of the sphere were not. In that geome- 
try the curvature is distributed inhomogeneuosly 
and therefore space-time symmetries present for 
the torus are violated. These effects have been 
studied and no consistent FSS at criticality has 
been observed. Instead the data result into an up- 
per bound on the discontinuity in the plaquette 
operator on the Wilson line, consistent with the 
earlier extrapolation in the vicinity of the TCP [4]. 

2. CQED on Almost Perfect Spheres 

The construction of an improved lattice action 
of CQED on a sphere is guided by the assump- 
tion, that certain non-universal features of the 
lattice discretization should become irrelevant in 
the vicinity of a field LGT fixed point. There- 
fore we may as well use the formulation of gauge 
theories on random lattices [5] and in close anal- 
ogy construct an improved action function, which 
should correct some of the defects of the above 



described naive discretization. In a first step we 
project sites, links and plaquettes of the mani- 
folds SH(L) and SH'(L) (dual) onto the surface 
of 4-dimensional spheres. We then consider an 
action of the form 
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J2^[(]cosiUp) + jcosi2Up)], 



(2) 



where each plaquette has now a weight, which is 
given by the ratio of its dual plaquette area di- 
vided by the plaquette area. A triangularization 
making use of euclidean distances approximates 
the plaquette areas Ap and A'p. Note that this 
lattice gauge theory behaves correctly in the naive 
(small coupling) continuum limit. In addition this 
mapping of the regular manifold SH(L) onto the 
sphere preserves the intrinsic lattice dimensional- 
ity Dh = 4. 

3. Results 

We have studied the above theory (2) at two 
values of the parameter j = and j = —0.2. 
We simulate the theory using a 3-hit Metropolis 
update (including overrelaxation steps at 7 = 0) 
on a sequence of lattice sizes constructed in the 
above mentioned way from SH(L) lattices with 
_L = 4, 6, 8 for both values of 7, and _L = 10 at 
7 = 0. In the vicinity of the finite volume phase 
transition we simulate at fixed value of 7 sev- 
eral /3-values and employ the Ferrenberg Swend- 
sen re-weighting method for the construction of 
continuous thermodynamic quantities as function 
of l3. For each lattice size we constructed typi- 
cally (0.5 — 1) X 10® configurations. We consider 
moments of the plaquette energy i.e., the specific 
heat C'v and other cumulants like the Binder cu- 
mulant VpcL and a quantity U4 (see e.g. [4] for 
definitions). These quantities develop extrema at 
criticality and their locations define pseudocriti- 
cal couplings. 

Our numerical analysis indicates a critical be- 
havior of the theory at the 7- value of 7 = —0.2. 
Even at 7 = the behavior on the given lattice 
sizes is consistent with a critical point. To il- 
lustrate these statement we display in fig. 2 the 
probability distribution of the plaquette energy 
P(E) on a _L = 8 lattice at 7 = —0.2 at the 
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peak position of the specific heat. The variable 
E denotes the /3— term of the action (2). The dis- 
tribution function is consistent with a Gaussian 
form, which has to be contrasted with the data 
of fig. 1. Such Gaussian forms are observed for 
all our lattice sizes at 7 = —0.2, small deviations 
from this form show up for the data at 7 = 0. In 
fig. 3 we display our data for the maxima of the 
specific heats in a double logarithmic plot as a 
function of the linear size scaling variable L. As- 
suming standard finite size scaling arguments at 
a critical point we expect Cv,max to scale cx . 
Assuming furthermore the hyperscaling relation 
a = 2 — Diy we expect a behavior cx L~~^ . The 
observed slopes are f = 1.47(20) and 1.18(30) for 
7 = and 7 = —0.2 respectively. These values 
are far away from the value 4 expected for a first 
order phase transition. We find i> = 0.36(2) for 
7 = and = 0.38(4) at 7 = -0.2. In addition 
we have checked the finite size scaling of pseudo- 
critical /3- values, which at criticality are supposed 
to scale l3pc(L) = (3pc{L = 00) -|- c/ L~ . The data 
are consistent with such scaling laws and values 
for v close to mentioned ones. 
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Figure 2. P{E) on sphere at 7 = —0.2. 

4. Conclusion 

Using a regularization of CQED on a sphere 
we observe a critical, continuous behavior of ther- 
modynamic quantities for the value of 7 = —0.2. 
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Figure 3. Maxima of the specific heat. 

The observed FSS behavior is consistent with the 
singular behavior of a critical point, i.e. it is 
suggested that the theory (2) exhibits a corre- 
lation length divergence there. At 7 = we ob- 
serve an "effective" critical behavior, but due to 
the limited lattice sizes we cannot draw a defi- 
nite conclusion. The measured i> values indicate 
a nonstandard situation, where the value v = ^ 
does not seem to be realized. This observation is 
connected to the question whether the underlying 
topological disorder transition constitutes a new 
universality class in AD LGT. We are currently 
accumulating more data on larger lattices and at 
7 = —0.5. Further analysis, including a study of 
zeros of the complex partition function, then will 
allow us to sharpen our statements. 
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